clc,clear
syms p h L a
b = a / (2*pi);
d = 25;
c0 = 1479/80*pi*2;
c1 = 1741/80*pi*2;
c2 = 1947/80*pi*2;
c3 = 2123/80*pi*2;
c4 = 2280/80*pi*2;
eqns = [
    0.5*h == (d+a)/2*(c2-c0) - b/2*(c2^2-c0^2)
    0.5*h == (d+a)/2*(c4-c2) - b/2*(c4^2-c2^2)];
vars = [h,a];
[h,a] = solve(eqns,vars)
